Wigner crystallization and metal-insulator transition of two-dimensional holes in 

GaAs at B=0 



Jongsoo Yoon, C.C. Li, D. Similar 1 , D.C. Tsui, and M. Shayegan 
Department of Electrical Engineering, Princeton University, Princeton, NJ08544 
1 Department of Condensed Matter Physics, Weizmann Institute, Rehobot 76100, Israel 

(July 14, 1998) 



00 

On 
ON 



We report the transport properties of a low disorder two-dimensional hole system (2DHS) in the 
GaAs/AlGaAs heterostructure, which has an unprecedentedly high peak mobility of 7 x 10 5 cm 2 /Vs, with 
hole density of 4.8 x 10 9 cm~ 2 <p<3.72xl0 10 cm~ 2 in the temperature range of 50mK<T<1.3K. From their 
T, p, and electric field dependences, we find that the metal-insulator transition in zero magnetic field in 
this exceptionally clean 2DHS occurs at r s = 35.1 ± 0.9, which is in good agreement with the critical r„ 
for Wigner crystallization r s c = 37 ± 5, predicted by Tanatar and Ceperley for an ideally clean 2D system. 

PACS numbers: 71.30.+h, 73.20.Mf, 73.20.Dx 
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It is well known that the ground state of a two- 
dimensional (2D) electron or hole system realized in semi- 
conductors at B=0 is expected to be a Wigner crystal 
(WC) in the low density limit when carrier-carrier in- 
teraction becomes dominant over the kinetic energy of 
the carriers. Tanatar and Ceperley |jj from their Monte 
Carlo simulation for ideally clean systems find that the 
critical density, expressed in terms of the dimension- 
less parameter r s , is r s = 37 ± 5. We note that r s 
is the Coulomb energy to Fermi energy ratio, given by 
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/4ttH ee Q , where m* is the effective 



mass, p the charge carrier density, and e the dielectric 
constant. Disorder, which is always present in any exper- 
imental system, also plays important roles. In the clean 
limit, it is expected to pin the WC, rendering the system 
into an insulator. In the highly disordered dirty limit, 
the Coulomb interaction is negligible, or can be treated 
as a small perturbation; the system is better described 
by single particle localization. In between, when both in- 
teraction and localization are equally important, neither 
the electronic processes nor the nature of the phase is 
known and there is no clear physical picture for the sys- 
tem at the present time. Also, it is believed that if the 
disorder is not too strong, it can enhance Wigner crys- 
tallization. Chui and Tanatar j|] find from their Monte 
Carlo studies that the disorder they put in their model 
can stabilize the WC phase at r s as low as 7.5. 

On the experimental side, a number of publications 
[ ^ ioJIoQiHeIJISI ' smce the discovery by Kravchenko et al. 
[3| of a 2D metal- insulator transition (MIT) in high mo- 
bility Si metal- oxide-semiconductor field-effect transis- 
tors (Si MOSFET's) at r s « 10, have reported transport 
studies in the large r s regime of 2D electron or hole sys- 
tems in several different semiconductor heterostructures. 
These studies generally do not address the properties in 
the insulating phase. They focus mostly on the nature of 
the metallic phase and the scaling behavior near the MIT. 
Their results agree with the conclusion of Kravchenko et 
al. pf that there is a MIT in high mobility samples and 



the critical density for the transition varies from r s 5 
in Sio.88Geo.12 to r s w 28 in p-GaAs. The origin and 
nature of the insulating phase, however, can be inferred 
from earlier Si-MOSFET experiments, where the low mo- 
bility samples show a cross-over from weak localization to 
strong localization as the carrier density is decreased. In 
high mobility samples, which show MIT's at r s w 10, Pu- 
dalov et al. |ll[] first studied the transport characteristics 
in the insulating phase and attributed them to those of a 
pinned WC. However, Shashkin et al. |12| and Mason et 
al. subsequently find that single particle localization 
gives a better description and the transport is via hop- 
ping in the presence of a Coulomb gap. It thus appears 
that in the regime where disorder and localization must 
be treated on equal footing, transport and I-V character- 
istics do not allow us to identify whether the insulating 
phase is a WC. A more desirable approach to the prob- 
lem has to come from clear identification of a 2D system 
in the clean limit, where the physics is expected to be 
simpler and theoretical predictions more reliable. 

In this letter, we report transport measurements on 
a 2D hole system (2DHS) in GaAs/Al x Gai_ x As, which 
has an unprecedentedly high peak mobility, /i p = 7 x 
10 5 cm 2 /Vs. This mobility, approximately twice time 
higher than previously studied 2DHS, allowed us to reach 
p as low as 4.8 x 10 9 cm -2 . From a systematic study of 
density, temperature and bias dependences of the hole 
transport, we find a MIT at r s = 35.1 ± 0.9, which is in 
good agreement with r s = 37 ± 5 predicted for the WC 
transition by Tanatar and Ceperley. We also examine the 
dependence of r s at the MIT, from all 2D electron and 
hole systems that have been reported to show MIT's, on 
the disorder in the sample and show that our sample is 
indeed in the clean limit. In the rest of this paper, we 
present this result and also describe in detail related new 
transport characteristics from our 2DHS. 

The 2DHS is created in a modulation doped GaAs/ALj 
Gai_a;As heterostructure grown on the (311)4 surface of 
GaAs substrate by molecular-beam epitaxy. A backgate, 
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used to change the density, is situated 250pm away from 
the 2DHS, and does not screen the carrier-carrier inter- 
action. Since the holes are produced by modulation dop- 
ing, the background neutralizing charges are immobile. 
In contrast, the 2D carriers in the insulated gate field- 
effect transistors used in other experiments HEHillfl are 
capacitivcly induced by applying a voltage to a metallic 
gate located close to the 2D systems. In such transistors, 
the background neutralizing charges are mobile and the 
carrier-carrier interaction is screened. The sample was a 
rectangle of 4.2 x 0.4mm 2 cut along the [233] direction, 
after the wafer was thinned to 250/im. Resistivity (p) 
was measured by a four probe method. To avoid sample 
heating, current was limited to below InA, and voltage 
was measured across two contacts separated by 0.7mm 
located near the center. An ac technique using a lock-in 
at a frequency <10Hz was used whenever possible. At 
low densities where capacitive coupling became signifi- 
cant, a dc method was employed. The sample was cooled 
in the mixing chamber of a dilution refrigerator. Three 
samples cut from the same wafer were studied, and all 
showed MIT's with the same characteristics. However, 
the data reported in this Letter were collected over three 
cool-downs on the same sample. The sample was warmed 
up to room temperature between cool-downs. 

Fig. 1 shows the T dependence of p, measured in 
the zero current limit, in the range 50mK<T<1.3K, for 
4.8 x 10 9 < p < 3.72 x 10 10 cm" 2 or 16.0 < r s < 44.4. 
For p>2x 10 10 cm -2 (the two bottom traces), the sign 
of dp/dT, which has been usually used to distinguish 
metallic phase from the insulating phase, is positive 
in the entire T range, showing metallic behavior. For 
p<7.2x 10 9 cm -2 (the four top traces), dp/ dT is always 
negative exhibiting the characteristic of an insulator. In 
the intermediate range of p, however, there is a peak in 
p resulting from dp/dT<0 at high T and dp/dT>0 at 
low T. This peak starts to appear at p=1.98x 10 10 cm~ 2 ; 
it becomes more pronounced and shifts to a lower tem- 
perature as p decreases. The maximum is at T«0.7K at 
p=1.98 x 10 10 cm- 2 , and «0.18K at p=9.0 x 10 9 c m - 2 . 
The presence of such a peak makes it ambiguous to tell 
whether the 2DHS is in the metallic phase or insulat- 
ing phase using the sign of dp/dT as the criterion. For 
example, at p=9.0 x 10 9 cm -2 , if the measurements are 
made only down to 0.3K, one would conclude from the 
fact dp/dT is negative that the 2DHS is in the insulat- 
ing phase, while the data below 0.18K indicate that the 
2DHS is in fact in the metallic phase. 

However, this ambiguity can be avoided if we de- 
fine the phase as metallic when the I-V is superlin- 
ear (d 2 I/dV 2 <0) and insulating when it is sublincar 
(d 2 I/dV 2 >0). Fig. 2 is the differential conductance 
(dl/dV) plotted against V at p=9.0 x 10 9 cm~ 2 , taken at 
several temperatures below and above the p maximum, 
as marked by the arrows in the inset. In sharp contrast 
to dp/dT, which shows a metallic behavior below and 



an insulating behavior above T=0.18K, the dl/dV vs. 
V taken at all four temperatures show superlinear I-V's, 
i.e., dl/dV decreases with increasing V. Such a superlin- 
ear I-V is observed down to p=8.1 x 10 9 cm~ 2 . When p is 
reduced below 7.3 x 10 9 cm~ 2 , where p shows monotonic 
increase with decreasing T, the I-V becomes sublinear, in- 
dicative of an insulating phase. It should also be pointed 
out that identifying the metallic phase and the insulat- 
ing phase based on the sign of d 2 I/dV 2 gives exactly the 
same result as that obtained using the sign of dp/dT at 
our lowest T. 

Using the sign of d 2 I/dV 2 , we locate the MIT at 
p=7.7±0.4x 10 9 cm~ 2 , or equivalently at r s = 35.1±0.9. 
At the MIT, p is approximately 35kSl at T-^0, and E F 
is calculated to be 0.58K. The critical r s for the MIT 
in our sample is approximately four times that in high 
mobility Si MOSFET's and 25% larger than that in pre- 
viously studied highest mobility 2DHS in GaAs/AlGaAs, 
of which p p k 4x 10 5 cm 2 /Vs is approximately half of that 
in our 2DHS. It is clear from this result that the r s , at 
which the 2D MIT is observed, depends strongly on the 
disorder in the 2D system. In the inset of Fig. 1, we plot 
r s c , the 2D MIT critical r s , from all samples reported in 
the literature as a function of the scattering rate of the 2D 
carriers, defined by 1/r = e/(m*p p ), which gives a simple 
measure of the disorder and also takes into account of the 
differences in effective mass, m*, of the charge carriers in 
different material systems. It is clear that for 1/r >lx 
10 11 sec" 1 , where all the high mobility Si MOSFET data 
lie, r s c is approximately constant. The data from the 
2D systems in GaAs/AlGaAs heterostructures, which are 
known to have low disorder and show fractional quantum 
Hall effect characteristics in high magnetic fields, are all 
in the l/r<l x 10 11 sec _1 regime. In this low scattering 
rate region, r s c shows a steep increase with decreasing 
1/r and extrapolates to r s c «40 as 1/r— >0. Our data, 
shown as the filled circle in the figure is in the clean limit 
and the observed r s c = 35.1 ± 0.9 is in good agreement 
with the r s = 37 ± 5 for the WC transition predicted by 
Tanatar and Ceperley. 

In the clean limit, a WC is expected to be pinned, giv- 
ing rise to a depinning threshold voltage, which is char- 
acterized by a sudden increase in current as the applied 
voltage exceeds the threshold. In Fig. 3, the dV/dl 
taken at r s =44.4 is plotted against V. At 50mK, the 
dV/dl is nearly constant until V reaches about 1.5mV, 
beyond which it drops rapidly by an order of magni- 
tude at V«3mV. This is suggestive of the depinning 
threshold expected for a pinned WC. In the WC con- 
duction model proposed by Chui 14 1, where the WC 
is weakly pinned by the potential of remote dopants, 
transport is expected to be mediated by the creation 
of dislocation pairs through quantum tunneling. The 
energy needed to create a dislocation pair leads to a 
conduction threshold voltage. The observed threshold 
voltage of 1.5mV, corresponding to a field of 2.1V/m, 
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is in reasonable agreement with Chui's model, accord- 
ing to which, the depinning threshold field estimated for 
our sample is E d = 0.09nja 2 e 2 /ed 3 = 2.5V/m. Here, 
rii = 8 x 10 n cm~ 2 is the remote doping impurity density, 
a = (7rp) -1 / 2 , and the doping impurity setback distance 
d=1200 A. However, there are caveats. The pinning gap 
deduced from E d , A p — Ti(eE d /m*a) 1 /' 2 «2.7K, is not 
consistent with the observation that no well defined de- 
pinning threshold feature is seen when the thermal energy 
is still much smaller than the gap. In Fig. 3, the data 
show no clear threshold at T as low as 145mK. Moreover, 
as shown in the inset of Fig. 3, the I-V taken at T as low 
as 50mK is linear within the experimental resolution at 
voltages below 1.5mV. This ohmic conduction below the 
depinning threshold voltage cannot be explained by pos- 
sible backgate leakage current, which is much less than 
lpA. 

Finally, we turn to the overall symmetry in the non- 
linear transport properties of our low disorder 2DHS as 
its density is varied across the MIT. Shahar et al. jl5| re- 
cently reported a new symmetry across the quantum Hall 
liquid to Hall insulator transition. They find that the TV 
curves in the quantum Hall liquid phase can be mapped 
into the I-V curves in the Hall insulator phase by simply 
interchanging I with V. This I-V inversion symmetry can 
be viewed as reflecting the charge-flux duality symme- 
try in the composite boson description of the fractional 
quantum Hall effect. More recently, Kravchenko et al. 
j| and Simonian et al. [jl6| reported a similar reflection 
symmetry across the 2D MIT, though its origin is still 
unclear at the present time. We have made a detailed 
comparison of the I-V characteristics in the metallic and 
insulating phases in our sample, and find that they are 
qualitatively different. Instead of comparing I-V with V-I 
in the two phases, we compare their derivatives which are 
more sensitive to changes in the curvature. In Fig. 4(a) 
and 4(b), dl/dV in the insulating phase and dV/dl in the 
metallic phase are shown against I and V, respectively. 
In the metallic phase, as V increases, dV/dl initially in- 
creases, but stops increasing at Vw 0.15mV. Above this 
characteristic voltage, dV/dl stays constant to I as high 
as lOOnA, the highest current used in our measurements. 
If we assume a reflection symmetry with respect to a 
straight line in I-V corresponding to p=35kr2/sqr. at the 
MIT, the characteristic voltage of 0.1 5mV in the metal- 
lic phase implies a characteristic current of 2.5nA in the 
insulating phase, above which dl/dV should remain con- 
stant. It is clear from Fig. 4(a) that such a characteristic 
current is absent in all traces of the dl/dV vs. I in the 
insulating phase. We conclude that the TV reflection 
symmetry across the 2D MIT is not observed in our data 
in the ranges of I and V shown in Fig. 4. 

To summarize, we presented transport properties of a 
2DHS that undergoes a MIT at r s ° = 35.1 ± 0.9, which 
is unambiguously determined using the sign of d 2 I/dV 2 
rather than dp/dT. The dependence of r s c on disorder, 



from all 2D electron and hole systems that have been re- 
ported to show MIT's, indicates that our 2DHS is in the 
clean limit, and the value of r s c in our sample is in good 
agreement with r s = 37 ± 5 expected for the WC tran- 
sition by Tanatar and Ceperley. In the insulating phase, 
a conduction threshold that is quantitatively consistent 
with the WC conduction model of Chui is found, but the 
pinning gap estimated from the observed threshold field 
is too large compared with our observations. We also 
find that the TV reflection symmetry across the MIT is 
absent in our low disorder 2DHS. 
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Fig. 1. T dependence of p, from the top curve, at 
p=0.48, 0.55, 0.64, 0.72, 0.90, 1.02, 1.27, 1.98, 2.72, and 
3.72 in units of 10 10 cm" 2 . The inset shows the de- 
pendence of r s c on 1/r, from various 2D systems; this 
work (filled circle), 2DHS in GaAs/AlGaAs (open circles. 
|,g§), 2DES AlAs/GaAlAs (filled triangle, §), 2DHS 
in Sio.8 8 Ge .i2/Si (filled square, (k]]), and 2DES in Si- 
MOSFET (open triangles |,@,§). m*/m e =0.37, 0.46, 
0.18, and 0.19 for p-GaAs |l7|], n-AlAs @, p-Sio.8 8 Ge .i2 
||, and n-Si-MOSFET. 

Fig. 2. dl/dV in the metallic phase near the MIT 
(p=9.0x 10 9 cm -2 ) at four different temperatures marked 
by the arrows in the inset. Qualitatively same nonlinear 
I-V characteristics are seen at all four temperatures. In 
contrast, dp/dT shows a metallic behavior below and an 
insulating behavior above T=0.18K. 

Fig. 3. dV/dl at p=4.8xl0 9 cm -2 is plotted against 
V at T=50mK (open circles), 145mK (filled circles), and 
284mK (triangles). The inset shows the I-V curve at 
T=50mK. 

Fig. 4. (a) dl/dV vs. I in the insulating phase and 
(b) dV/dl vs. V in the metallic phase at 50mK. The 
numbers shown next to the symbols are the densities in 
units of 10 10 cm~ 2 . The corresponding rs values are in 
parenthesis. 
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